Dynamical scaling analysis of the optical Hall conductivity in the quantum Hall regime 
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Dynamical scaling analysis is theoretically performed for the ac (optical) Hall conductivity 
o'xyisp,^) as a function of Fermi energy ep and frequency to for the two-dimensional electron 
gas and for graphene. In both systems, results based on exact diagonalization show that axy{eF,ijj) 
displays a well-defined dynamical scaling, for which the dynamical critical exponent as well as the 
localization exponent are fitted and plugged in. A crossover from the dc-like bahavior to the ac 
regime is identified. The dynamical scaling analysis has enabled us to quantify the plateau in the 
ac Hall conductivity previously obtained, and to predict that the plateaux structure in ac is robust 
enough to be observed in the THz regime. 



Introduction — 

Dynamics of electrons in the integer quantum Hall ef- 
fect (QHE) is an interesting, hitherto not fully explored 
problem. Theoretically, the question is how the static 
Hall conductivity, which may be regarded as a topo- 
logical quantitjiii^, evolves into the optical Hall con- 
ductivity, especially in the THz regime where the rel- 
evant energy scale is the cyclotron energyi^"— Two of 
the present authors and Hatsugai have recently shown 
that the plateau structure in axyioj) is retained in the 
ac ('^ THz) regime in both the ordinary two-dimensional 
electron gas (2DEG) and in graphene (described as the 
massless Dirac model) , although the plateau height devi- 
ates from the quantized values in acj^ The numerical re- 
sult indicates that the plateau structure remains remark- 
ably robust against disorder, which can be attributed to 
an effect of localization which dominates the physics of 
electrons around the centers of Landau levels in disor- 
dered QHE systems. However, what is physically signif- 
icant is not a result for a specific sample size, but the 
scaling behavior, especially when the localization is rel- 
evant in disordered systems. For ac responses, we have 
to look into the dynamical scaling. Scaling analysis of 
localization-delocalization transition in the 2DEG QHE 
has been done for both the static longitudinal conduc- 
tivity axxisp), where £f is the Fermi energy)^ and for 
dynamical scaling properties of the longitudinal conduc- 
tivity axx{£F,'^)r' but the dynamical scaling for the Hall 
conductivity axy{sF-,uj) has not been properly addressed 
for both ordinary and graphene QHE. 

With this motivation, here we elucidate the dynami- 
cal scaling behavior of the ac Hall conductivity around 
the plateau to plateau transition to gain a deeper un- 
derstanding of the optical Hall effect and its robust step 
structures in the ac region. Namely, when we perform 
a scaling analysis for the plateau to plateau transition 
width W, the quantity depends on w. Physically, a new 
length scale, ~ i^^^, emerges at finite frequencies, 
where z is the dynamical critical exponent. We have 
performed the dynamical scaling analysis for both 2DEG 
and graphene QHE. The quantum Hall effect in graphene 
is unique in that a zero-energy Landau level (LL) exists. 



which has no counterpart in the QHE in 2DEG^. Thus 
the dynamic scaling is of special interest for the n = 
LL in graphene. 

Experimentally, scaling properties of crxxi^Fi^^) was 
investigated from w = 0^ up to the GHz regime^. Recent 
advances in optical measurements (e.g., Faraday rota- 
tion in magnetic fields) in the THz region have made the 
study of dynamical response functions feasibleiiii^. For 
graphene, optical properties begin to be studied, among 
which are experimental transmission spectra^^, or the- 
oretical examination of the cyclotron emission^. Thus, 
the physics of dynamical scaling in graphene QHE should 
be interesting in the THz regime. 

Here we shall show that: (i) The ac Hall conductivity 
obeys a well-defined dynamical scaling, (ii) There is a 
crossover in the scaling behavior from a dc-like regime 
to an ac regime, in the latter of which dominates the 
scaling. In the former L/^ dominates the scaling (where ^ 
is the localization length), while in the latter Li^/£, does, 
(iii) The dynamical critical exponent is found to be z ~ 2 
in both the 2DEG and graphene QHE systems as far as 
the potential disorder is concerned, (iv) The analysis 
enables us to estimate the plateau to plateau transition 
width W in the ac regime with Luj < L to assert that 
the Hall conductivity maintains the plateau structure at 
frequencies as high as a; ~ O.lwc, which, for a magnetic 
field of a few Tesla, covers the THz region. This is an 
experimentally testable statement. 

Formalism — 

For the ordinary QHE system as typically realized in 
GaAs/AlGaAs, the kinetic part of the Hamiltonian is 
Hq ~ ■^^{p + eA)'^, where m* is the effective mass of the 
electron, p = {px,Py) the momentum, and A the vector 
potential. Disorder is introduced by a random potential 
V{r) composed of Gaussian scattering centers of range d 
and density nimp placed on randomly chosen points Rj: 

l/(r) = ^w,exp(-|r-R,f /2d2)/(27rd2). (1) 
j 

For Uj we assumed a bimodal distribution Uj = ±u 
with random signs so that the broadened Landau level 
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is symmetric. A measure of disorder, i.e., the Landau 
level broadeningi^, is T = 2M[ni„ip/27r(^2 _^ (p)]^/^. Here 

we take d = 0.7£, where £ = y^fi/eB is the magnetic 
length, but the result does not change significantly for 
other choices of d. Diagonalization of the Hamiltonian is 
done for the subspace spanned by the five lowest LL's for 
L X L systems with L/£ varied over 25,30,35,40. With 
wave functions and energy eigenvalues ea at hand, the 
optical Hall conductivity is evaluated from the Kubo for- 
mula, 



i ^ ^ Cfc - Ca 



■ab-ba 
J X Jy 



■ab-ba 
Jy J X 



Eb- £a~ Tt-UJ £6 - Ea + Tt-UJ 



(2) 
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where j^^ is the current matrix element^. 

The conductance is then averaged over a few thou- 
sands samples with different disorder potential realiza- 
tions. The averaged conductivity is hereafter denoted by 
the same symbol (Jxyi^F, w). For the scaling analysis the 
calculation done for varied sample size L, energy £p and 
frequency w. 

For graphene QHE, we employ the two-dimensional ef- 
fective Dirac model. 



H = vpcr ■ TT + V{y), 



(3) 



where cr = {ax,(Jy) is the Pauli matrices, tt = p -I- eA, 
and V{y) the random potentiali^. The selection rule for 
the current matrix elements in the Dirac model (|n| ^ 
|n| ± 1 with n the Landau index) is distinct from that 
(n o n ± 1) for 2DEG. 

Raw results for the Optical Hall Conductivity — 
The optical Hall conductivity Uxyi^F,^) as a func- 
tion of the Fermi energy ep and frequency u is dis- 
played for the 2DEG (FiglUa)) and graphene (FiglUb)) 
QHE systems. The density of states (DOS) for Lan- 
dau levels (insets of Fig[lja,b)) confirms that the Lan- 
dau level broadening is ~ F. For each value of ef, 
the frequency-dependence is the Hall conductivity can 
be recognized as the cyclotron resonance. The 2DEG 
has one resonance at cyclotron frequency uic (FigUJa)), 
while the graphene QHE system exhibits a series of reso- 
nances, which correspond to the Dirac QHE selection rule 
\n\ o |n| ± 1 (FiglTJb)) for the non-uniform set of Lan- 
dau levels (cx \/n). Away from a resonance, a step-like 
structure is seen in Uxy [ef,^) as a function of e^? even for 
finite values of w, although the step heights are no longer 
quantized. We can attribute this behavior of crxy{eF:^) 
to the localization property of electrons in QHE systems: 
As long as w <C Wc,r, the nature of the mobility gap is 
maintained and Uxy remains flat. If we more closely in- 
sepct the width W of the plateau to plateau transition, 
W for a finite w in Figlljd) is seen to be greater than in 
the static case in Figlljc), although still narrower than 

r. 



FIG. 1: axy{eF,ijj) plotted against Fermi energy ef and fre- 
quency uj for (a) 2DEG, and (b) graphene quantum Hall sys- 
tems for a pontential disorder with P = OAhujc- Insets are 
density of states for the lowest (n = 1) LL (a) and n = LL 
(b). (c)The static Hall conductivity (Txy{eF) plotted against 
ef for the graphene QHE system of sizes L = 25,40. (d) 
Optical Hall conductivity Gxy{eF plotted against ef for 
the graphene QHE system of sizes L = 25, 40 for a fixed 
UJ — Qujc/L^)- The solid lines in (c) and (d) represent fit- 
ting with eqn.(l5l). 



We now move on to the sample size dependence of the 
widths W{e,Ljj,L) for the static and dynamic Hall con- 
ductivities. For the static Hall conductivity CTa;y(e, 0,L), 
we confirm the standard picture, where the plateau to 
plateau transition width becomes narrower with the sam- 
ple size L as seen in FigHJ In the thermodynamic limit, 
almost all the wave functions are localized, where the lo- 
calization length diverges like ^ ^ 1/|£f — ed" toward the 
center of the LL at e = Ect^ For finite systems the states 
whose localization length ^ is larger than the system size 
L are effectively extended, and contribute to the longitu- 
dinal conductivity and the plateau to plateau transition. 
This suggests the behavior W ^ L"^/"^. 

Scaling Analysis — We are now in position to look at 
the dynamical scaling analysis of the optical Hall conduc- 
tivity (Jxyi^F, ^, L) and the width W{u)^ L) of the plateau 
to plateau transition. We expect that the W increasing 
with oj and decreasing with L may be captured with some 
scaling function. For that we have to quantify the width 
W , or the steepness (cx l/H^) of the transition by fitting 
<yxy{£F,^, L) around the transition region for a given LL 
to some function of e^? for each value of a;. To describe 
the transition Uxyli-e^jh) = ^ 1 in the 2DEG QHE 
we take 



Gxy{eF-,<^,L)2T)^G = ^ + ^tanh 



Ef - ^hiVc 



W{lu,L) 



(4) 



while for the transition axy/{—e'^/h) = —1 — !■ 1 in 
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graphene QHE we take 



0'a:j,(eF, L)gi.aphGno — tailh 



SF 



WiLo,L) 



(5) 



The quality of fitting of the plateau to plateau transition 
by the tanh function is quite satisfactory^! as can be seen 
in Figlli;c,d). 

Dynamical scaling analysis for axy{eF,^, L) is carried 
out in a similar manner as that for the longitudinal 
conductivityi^. In this ansatz the optical Hall conductiv- 
ity is regarded to depend on Fermi energy and frequency 
only through the ratios L/^ and L^/£^. Here we have the 
localization length, ^ ^ l/\eF — where £c is the criti- 
cal energy which coincides with the center of the LL, and 
Then the dynamical scaling ansatz for the 
optical Hall conductivity reads 



axy{eF,uj,L) = -F{{eF ~ ec)L^/'',ujL'), (6) 

where is a universal scaling function. This implies that 
the width of the plateau to plateau transition scales as 



(7) 



where / is a universal function deduced from F. The 
first factor on the right-hand side makes the plateau to 
plateau transition width narrower for larger systems and 
dictates the dc scaling, while the second factor / de- 
scribes the dynamical scaling. 

In Figl2]we show the scaling of inverse transition width, 
l/W{uj,L), for the 2DEG QHE system. By examin- 
ing first the inverse width 1/W{uj = 0) for the static 
case against system size L in Figl^^b), we obtain the 
localization critical exponent ly from logl/W(a; = 0) = 
l/iy\ogL+f{0), with the resuU v = 2.1±0.2. This agrees 
with the accepted value of the static critical exponent in 
the integer QHE, albeit slightly smaller. 

On the other hand, the frequency dependence of the in- 
verse width for a fixed system size L in (Figl^Ja)) clearly 
exhibits that there are two regions: In the first region 
l/W stays nearly constant up to some critical frequency 
that depends on the system size L, while in the other the 
quanitity begins to decrease monotonically with lj. In 
the latter region, l/W assumes similar values for all the 
sample sizes studied here as shown in (Figl^Ja)). We can 
indeed notice that, in the first region we have L < Li^, 
while in the second L > L^i- If we inspect eqnl?! and as- 
sume a power-law form for the scaling function /, we can 
see that the inverse width in the second region should 



take a (L-independent) form, 1/W{uj) 



OC UJ 



-1/^". Cal- 



culation of the dynamical exponent z should be done 
for the critical region (i.e. for the transition width not 
too large), so that we do this around the crossing re- 
gion where l/W begins to decrease. This happens, typ- 
ically, for UJ < 0.002a;c. With a least square fitting 
of logl/W {ui) ~ const — ^logoj, the dynamical criti- 
cal exponent for the 2DEG QHE system is obtained as 
z= 1.8 ±0.2 . 
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FIG. 2: Dynamical scaling analysis for the 2DEG QHE sys- 
tem with r = OAhujc and L — 25 — 40: (a) The inverse width 
l/W plotted against the frequency u, (b) the inverse width 
l/W plotted against the sample size L, from which the local- 
ization exponent = 2.1 ± 0.2 is obtained, and (c) rescaled 
inverse width -^L~^^'^ plotted against the rescaled frequency 
ujL" with a fitted dynamical critical exponent z = 1.8 ± 0.2. 
W,uj,L are measured, respectively, in units of huJc,ujc,£- 



If we now turn to the graphene QHE system, the same 
analysis is performed based on FiglH The frequency 
dependence of the transition width is rather similar to 
that for the 2DEG system, as far as the potential dis- 
order assumed here is concerned. The localization ex- 
ponent v and and the dynamical critical exponent z for 
the graphene system are determined as v = 2.1±0.1 
z = 1.8 ± 0.2, which coincide, within numerical errors, 
with those for the conventional QHE system. This sug- 
gests that the two systems are in the same universality 
class. As far as the dynamical exponent is concerned, 
it has been argued by Hikami and Wegner— that, when 
the density of states is an analytic function of energy 
at a critical point, then z = d [d: spatial dimension, 
which is 2 in the present case). Thus the fact that we 
get z w 2 for the 2DEG QHE is commensurate with this 
conjecture. On the other hand, the density of states for 
Dirac fermions {p{E) ^ \E\ for the clean system) is non- 
analytic, for which one might expect a different behavior 
in graphene. The absence of deviation in the value of z in 
Dirac fermions here should come from the fact that the 
presence of disorder smears the Dirac cone structure in 
the density of states to make it smooth (FiglTJb), inset). 
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FIG. 3: Dynamical scaling analysis for the graphene QHE 
system with F = OAhijJc and L = 25 — 40: (a) The inverse 
width 1/W plotted against the frequency lo, (b) the inverse 
width 1/W plotted against the sample size L, from which 
the localization exponent = 2.1 ± 0.1 is obtained, and (c) 
rescaled inverse width -^L"^^'^ plotted against the rescaled 
frequency ljL'' with a fitted dynamical critical exponent z — 
1.8 ±0.2. 



Having derived the static and dynamic exponents, we 
can now actually plot the scaling of the rescaled inverse 
width against rescaled frequency. This is displayed in 
Figlllc) for 2DEG, and FigElJc) (for graphene QHE). It 
can be judged that the scaling fit is quite good, which 
indicates that the form of the universal function assumed 
in eqn(I7]) is adequate. In the scaling plot we can see 
more clearly the first region with a constant ■^L~^^'^ for 
small LuL^, and the second region with a monotonously 
decreasing 1/W for larger wL^. 

Intuitively, we can elaborate as follows: The dynami- 
cal response of the QHE system is governed by the mag- 
nitude of the localization length relative to two length 
scales: the system size L, and the length which is 



the distance over which an electron travels during one 
cycle, l/w, of the ac field. Since the localization length 
diverges as ^ ^ |e — ed^ near the center of LL, and the 
states contributing to (Jxy{uj) should be those that simul- 
taneously satisfy £, > L and ^ > L^^ , the transition width 
is determined by the smaller length scale L or Luj. In the 
static limit w = 0, with — > oo, the system size L de- 
termines the transition width W—. When uj is increased, 
decreases. In the low enough frequency region, one 
still has > L so that the transition width continues 
to be determined by the system size. When < L for 
higher frequencies, however, W begins to be governed 
by Li^ , and the transition width broadens monotonically 
with frequency. For even higher frequencies, L^^ becomes 
so small that the system is far from the critical region, 
and departure from the scaling should occur. 

From the functional form of f{ujL^') in eqn©, the fre- 
quency for which uL^ ^ 1, corresponds to the crossover 
region where the two regions overlap, that is, Li^ L, 
or ^ 1/ijj^l" . With the dynamical scaling argument 
with z ~ 2, we end up with ~ l/ix)^!"^ ~ t^l'^ where 
t is the diffusion time. Since square-root time evolution 
is a characteristic of diffusion processes, the dynamical 
response behavior indicates that the present disordered 
system is diffusive. 

The dynamical scaling here enables us to give an esti- 
mate of the transition width W in the THz region (with 
typically uj ~ O.lwc)- In the L^j-dominated regime, one 
obtains W (x w^/^''. The proportionality constant can 
be read out from the numerical result. Fig. [21[c), so that 
W/hujc - 0.2a;i/^'' ~ 0.1 at w = O.lwc. This implies 
that the plateau structure remains robust up to the THz 
region, so that experimental measurements should be fea- 
sible. 

To summarize, dynamical scaling analysis of the op- 
tical Hall conductivity in 2DEG and graphene quantum 
Hall systems have been carried out, while previous stud- 
ies have focused mainly on the longitudinal conductivity. 
The dynamical critical exponent z ~ 2 implies that the 
system is in the diffusive limit. The dynamical critical 
exponent z is found to be similar between the 2DEG and 
the graphene QHE system, but we have to re-emphasize 
that this is as far as the pontential disorder taken here is 
concernded. It is well-known^ii that the preservation or 
otherwise of the chiral symmetry in disordered graphene 
has a profound effect on the n = graphene Landau 
level. Thus it is an interesting future problem to look 
into this effect in terms of the dynamical scaling. 
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